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PREFACE 


The  study  described  herein  was  sponsored  by  the  Office,  Chief  of 
Engineers,  U.  S.  Army  (OCE),  as  a  part  of  the  Mobility  and  Weapons 
Effects  Technology  RDT&E  Project  No.  UAT62719AT40,  Work  Unit  001,  "Air¬ 
field  Pavement  Design  and  Parametric  Sensitivity  Analysis,"  and  Work 
Unit  003,  "Rigid  Airfield  Pavement  Load-Deformation  Response  Analysis." 

This  report  is  Report  3  of  a  three-report  series  concerning  the 
computer  programs  WESLIQID  and  WESLAYER,  which  provide  for  analysis  of 
rigid  multicomponent  pavements  with  discontinuities  on  liquid  founda¬ 
tions  (WESLIQID)  and  on  linear  layered  elastic  solids  (WESLAYER).  This 
report  is  a  user's  manual  for  WESLAYER.  Report  1  presented  the  theo¬ 
retical  background  of  the  programs  and  numerical  results  and  discussed 
the  capability  and  logic  of  the  two  programs.  Report  2  was  a  user's 
manual  for  the  WESLIQID  program. 

The  study  was  conducted  by  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  Geotechnical  Laboratory  (GL),  under  the  gen¬ 
eral  supervision  of  Dr.  Don  C.  Banks,  Acting  Chief,  GL;  Dr.  Paul  F. 
Hadala,  Assistant  Chief,  GL;  and  Mr.  Alfred  H.  Joseph,  Chief,  Pavement 
Systems  Division  (PSD),  GL.  Dr.  Yu  T.  Chou,  PSD,  was  in  charge  of  the 
study  and  is  the  author  of  the  report.  Professor  Y.  H.  Huang  of  the 
University  of  Kentucky,  who  originally  developed  the  computer  programs, 
assisted  in  the  study. 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Com¬ 
manders  and  Directors  of  the  WES  during  this  study  and  the  preparation 
of  this  report.  Mr.  Fred  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


Mult iply 

..  B y. 

To  Obtain 

Fahrenheit  degrees 

0.555 

Celsius  degrees  or  Kelvins* 

feet 

0.301*8 

metres 

inches 

2.5*+ 

centimetres 

pounds 

(force) 

1*  .UU8222 

newtons 

pounds 

inch 

(force)  per 

cubic 

0.2711* 

megapascals  per  metre 

pounds 

(force)  per 

inch 

175.1268 

newtons  per  metre 

pounds 

inch 

(force)  per 

square 

6.891*757 

kilopascals 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  =  0.555  (F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  =  0.555  (F  -  32)  +  273.15. 
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STRUCTURAL  ANALYSIS  COMPUTER  PROGRAMS  FOR  RIGID  MULTICOMPONENT 
PAVEMENT  STRUCTURES  WITH  DISCONTINUITIES— 

WESLIQID  AND  WESLAYER 


MANUAL  FOR  THE  WESLAYER  FINITE  ELEMENT  PROGRAM 


PART  I :  INTRODUCTION 

Background. 

1.  The  U.  S.  Army  Corps  of  Engineers  (CE)  has  realized  for  many 
years  that  much  of  the  maintenance  of  rigid  pavements  is  associated  with 
cracks  and  joints.  Current  CE  rigid  pavement  design  procedures  have 
certain  limitations  that  were  imposed  by  the  state  of  the  art  at  the 
particular  stage  of  development.  During  the  development  of  the  proce¬ 
dure,  it  was  necessary  to  make  simplifying  assumptions  and,  in  many  in¬ 
stances,  to  ignore  the  effects  of  cracks  and  joints.  Since  the  advent 
of  high-speed  computers  and  the  development  of  the  finite  element  method, 
a  more  comprehensive  investigation  of  the  state  of  stress  at  pavement 
joints,  cracks,  and  other  locations  in  multicomponent  pavement  structures 
than  previously  possible  can  now  be  achieved.  Consequently,  a  better 

and  more  reasonable  design  procedure  may  be  developed  for  rigid  pavements. 

Purpose 

2.  The  development  of  the  finite  element  programs  and  the 
analysis  of  computed  results  are  presented  in  Report  1  of  this  series. 

This  report  presents  a  user's  manual  for  a  computer  program  named 
WESLAYER.  The  program  computes  the  state  of  stress  in  a  rigid  pavement 
supported  on  a  layered  elastic  foundation,  as  well  as  in  the  supporting 
subgrade  soils. 


Scope 

3.  The  computer  program  is  described  in  this  report  to  give 
users  a  concise  picture  of  the  program  without  reference  to  Report  1. 
The  logic  of  the  programming  is  explained  with  the  illustration  of  a 
flowchart.  An  input  guide  to  the  computer  program  is  given,  and  four 
example  problems  are  presented  to  illustrate  input  procedures  for  use 
of  the  computer  program.  The  outputs  of  the  example  problems  are 
explained. 
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PART  II:  PROGRAM  DESCRIPTION 


k.  This  report  describes  a  finite  element  computer  program 
named  WESLAYER  programmed  for  the  analysis  of  concrete  pavements  sub¬ 
jected  to  multiple-wheel  loads.  The  program  is  developed  for  a  subgrade 
soil  represented  as  a  linear  layered  elastic  solid.  Any  number  of  layers 
can  be  accommodated.  Because  of  the  assumption  of  an  elastic  soil,  ver¬ 
tical  force  at  one  nodal  point  in  the  subgrade  causes  vertical  movements 
at  all  other  nodes  and  vice  versa.  This  behavior  is  different  from  that 
of  a  liquid  foundation  (WESLIQID).  In  the  liquid  foundation,  the  verti¬ 
cal  force  or  deflection  at  one  nodal  point  in  the  subgrade  is  not  af¬ 
fected  by  those  at  other  nodes . 

5.  The  program  determines  stresses  and  displacements  in  the 
pavement  and  in  the  supporting  subgrade  soil  due  to  loads  and  tempera¬ 
ture  warping.  Part  of  the  pavement  can  be  out  of  contact  with  the  sup¬ 
porting  subgrade  before  applying  the  load  and  the  temperature  gradient, 
and  the  program  determines  the  condition  of  contact  at  each  nodal  point 
after  the  application  of  loads  and  a  temperature  gradient.  Input  data 
of  the  programs  include  (a)  the  property  and  geometry  of  the  pavement 
and  subgrade  soil,  (b)  the  magnitude  and  distribution  of  the  loads,  (c) 
the  temperature  gradient,  (d)  gaps  under  the  pavement  at  certain  nodal 
points,  if  any,  and  (e)  joint  and  crack  conditions. 

6.  Multiple-wheel  loads  can  be  input  and  the  number  of  wheels 
is  not  limited.  Because  of  the  large  computer  storage  space  required, 
the  program  can  handle  only  two  slabs,  except  for  the  special  option 
where  a  four-slab  pavement  system  is  loaded  symmetrically  at  the 
pavement's  center.  At  the  joint,  the  program  considers  only  the  shear 
transfer  and  assumes  the  moment  transfer  to  be  zero. 
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PART  III:  PROGRAM  APPROACH 


7.  The  storage  space  required  for  the  program  depends  on  the 
total  number  of  elements  used  in  the  problem.  An  iteration  scheme  is 
used  in  the  program  so  that  the  computation  is  made  only  for  one  slab 
at  each  time.  This  scheme  results  in  a  great  reduction  in  computer 
time  because  the  number  of  equations  to  be  solved  each  time  is  reduced 
to  only  one  slab.  Two  series  of  iterations  are  involved  in  the  program: 
one  is  with  respect  to  subgrade  contact  and  the  other  is  with  respect 

to  load  transfer  across  the  joint. 

8.  In  the  iteration  with  respect  to  subgrade  contact,  the  con¬ 
tact  condition  at  each  node,  i.e.,  whether  the  slab  and  subgrade  are  in 
contact  or  not,  is  first  assumed,  and  the  iteration  with  respect  to  load 
transfer  proceeds  until  either  the  convergence  criteria  (DEL  in  Item  11 
of  the  input  guide.  Table  1*)  are  satisfied  or  the  maximum  allowable 
number  of  iterations  ( ICL  in  Item  10  of  Table  l)  is  reached.  At  this 
stage,  the  resulting  contact  condition  is  determined.  If  some  nodes 
originally  assumed  in  contact  are  found  out  of  contact  or  some  nodes  as¬ 
sumed  out  of  contact  are  found  in  contact,  the  newly  found  contact  con¬ 
dition  is  assumed,  and  the  process  is  repeated  until  the  same  contact 
condition  is  obtained.  This  can  usually  be  achieved  within  only  a  few 
iterations.  The  only  control  by  the  user  is  to  specify  the  maximum 
number  of  iteration  cycles,  NCYCLE  .  If  NCYCLE  =  1  ,  the  contact 
condition  between  the  slab  and  subgrade  is  known  a  priori ,  and  no 
iterations  are  needed. 

P.  T n  the  iteration  with  respect  to  load  transfer  across  the 

joint,  tt  ?  computation  is  made  successively  between  the  left  and  right 
slabs.  The  vertical  deflections  are  used  for  checking  convergence. 

10.  The  program  carries  out  the  computations  (shown  in  the 
flowchart.  Figure  l)  in  the  following  sequence: 


*  The  input  guide  (Table  1)  appears  in  Part  IV,  where  it  is  discussed 
in  detail. 
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Figure  1.  Flowchart  for  WESLAYER  program 
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a_.  Generate  stiffness  matrix  for  each  element  (of  the 
slabs)  and  then  superimpose  them  to  form  an  overall 
stiffness  matrix. 

b.  Compute  the  element  flexibility  matrix  of  the  subgrade 
and  then  superimpose  them  to  the  matrix  of  the  slab  to 
form  the  overall  flexibility  matrix. 

£.  Invert  the  flexibility  matrix,  using  the  inversion 
subroutine  SMIN  ,  to  obtain  the  stiffness  matrix  of 
the  subgrade. 

d_.  The  stiffness  matrix  of  the  subgrade  is  added  to  the 
stiffness  matrix  of  the  slab  to  cover  the  stiffness 
matrix  of  the  pavement  system. 

e_.  Store  the  stiffness  matrix  adjacent  to  joints  for  later 
use. 

_f.  If  it  is  known  that  gaps  exist  under  certain  nodes  in 
the  subgrade  soil,  the  gaps  are  read  into  the  program 
to  combine  them  with  the  computed  curls  of  the  slabs 
due  to  temperature  warping  to  form  the  initial  subgrade 
contact  condition. 

£.  Determine  the  nodal  reactive  condition  based  on  the 
subgrade  contact  condition. 

h.  If  externally  applied  loads  are  considered,  the  uni¬ 
formly  applied  loads  are  distributed  to  the  adjacent 
nodes  using  statics. 

_i.  Compute  the  displacements  of  the  left  slab,  assuming 
that  there  is  no  shear  and  moment  transfer  along  the 
joints,  i.e.,  that  slab  1  has  four  free  edges. 

j_.  Impose  deflections  along  the  joint  from  the  left  slab 
to  the  right  slab  and  compute  the  displacements  of  the 
right  slab.  This  is  done  with  a  fixed  boundary  condi¬ 
tion  at  the  joint. 

k.  Determine  the  nodal  forces  in  the  left  slab  due  to  the 
deflections  of  the  subgrade  at  each  node  in  the  two 
slabs.  (Note  that  this  process  is  not  needed  in  the 
WESLIQID  program  in  that  the  reactive  force  at  a  certain 
node  depends  only  on  the  deflection  at  that  node  but 
not  elsewhere . ) 

1_.  Compute  the  vertical  nodal  forces  along  the  joint  in 
the  left  slab  due  to  the  deflection  of  the  right  slab. 

m.  Compute  the  deflections  of  the  left  slab  again.  The 
values  are  compared  with  those  computed  in  the  previous 
iterative  cycle.  If  the  convergence  criteria  are  not 
met,  the  nodal  forces  in  the  right  slab  due  to  the  de¬ 
flections  of  the  subgrade  at  each  node  of  both  slabs  are 
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computed  and  the  process  between  steps  J  and  m  are 
repeated  until  the  deflections  at  the  preselected  nodes 
converge  to  a  specified  tolerance. 

n_.  Once  a  convergent  solution  has  been  obtained  or  the  maxi¬ 
mum  allowable  number  of  iterative  cycles  has  been  reached 
( ICL  in  Item  10  of  Table  l),  the  signs  of  the  deflec¬ 
tions  at  each  node  are  compared  with  those  of  the  ini¬ 
tial  (or  the  previous)  subgrade  contact  condition.  A 
change  of  sign  at  any  node  indicates  that  the  contact 
condition  at  these  nodes  has  changed.  Based  on  the 
renewed  subgrade  contact  condition,  the  computational 
process  from  steps  a  to  m  is  repeated.  The  iteration 
process  stops  when  either  the  contact  condition  ceases 
to  change  or  the  maximum  allowable  number  of  iterations 
(NCYCLE  at  Item  U  of  Table  l)  has  been  reached.  Note 
that  when  the  subgrade  contact  condition  is  changed,  the 
subgrade  stiffness  is  also  changed. 

£.  Once  the  subgrade  contact  condition  is  unchanged,  the 

computational  process  from  steps  j  through  m  is  repeated 
once  more  with  a  refined  convergence  criterion.  The 
controlling  variables  in  the  program  are  ICLF  in 
Item  10  and  DELF  in  Item  11  of  Table  1 . 

£.  The  stresses  at  selected  nodal  points  are  computed 
based  on  the  curvature  of  the  deflected  slab,  i.e., 
the  nodal  displacements. 

£.  Compute  stresses  and  deflections  in  the  subgrade  if  so 
desi red . 

_r .  Note  for  a  single  slab,  i.e.,  NSLAB  =  1  ,  steps  j  to  m 
are  skipped  if  a  full  bandwidth  NB  is  used. 

11.  For  a  single  slab,  steps  j  and  1  are  neglected;  i.e.,  the 
left  slab  is  the  only  slab. 

12.  In  computing  the  subgrade  reactive  forces  at  each  nodal 
point,  Boussinesq's  homogeneous  analysis  and  the  layered  elastic  theory 
are  used  to  formulate  the  flexibility  matrix  for  single-layer  and 
multilayer  subgrade  soil,  respectively.  In  the  multilayer  case,  the 
layered  elastic  theory  is  not  used  directly  at  every  nodal  point  in 
order  to  save  computer  time;  rather,  the  theory  is  used  to  compute  the 
deflections  at  21  different  offset  points,  and  interpolation  subroutines 
are  used  to  interpolate  the  deflections  at  node  j  due  to  the  load  at 
node  i . 
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PART  IV:  OPERATION  OF  THE  PROGRAM 

General  Discussion 

13.  As  with  other  numerical  procedures  for  solving  structural 
problems,  the  accuracy  of  the  finite  element  method  depends  greatly  on 
the  correct  use  of  the  technique.  While  the  computational  cost  and 
storage  space  increase  drastically  with  an  increased  number  of  elements, 
the  program  does  require  a  reasonable  number  of  elements.  The  element 
size  should  be  smaller  near  the  loads  and  joints  where  stresses  are 
transferred  to  another  slab.  In  some  cases,  the  minimum  number  of 
elements  for  a  particular  problem  has  to  be  determined  by  a  trial-and- 
error  procedure.  It  was  found  that  an  insufficient  number  of  elements 
can  cause  the  solution  to  diverge;  this  is  particularly  true  when 
temperature  warping  is  considered  and  gaps  exist  under  the  pavement. 
Also,  users  should  be  aware  that  the  aspect  ratio  of  an  element,  defined 
as  the  ratio  of  the  larger  dimension  to  the  smaller  dimension  of  a 
rectangular  element,  should  not  exceed  four  or  five  to  one.  It  is 
always  good  practice  for  the  beginning  user  of  this  program  to  familiar¬ 
ize  himself  with  the  program  by  using  different  numbers  of  elements  for 
a  particular  problem  and  then  comparing  the  results. 

1*4 .  The  input  guide  for  the  program  is  presented  later  in  this 
Part.  Special  features  in  the  correct  and  efficient  use  of  the  program 
are  presented  and  discussed  in  the  following  paragraphs . 

Dimension  requirements 

15.  The  method  developed  in  this  program  can  be  applied  only  to 
two  slabs.  The  dimensions  of  C  *  and  G  vary  with  the  number  of  ele¬ 
ments  and  the  half  bandwith  and  can  be  computed  as 

NSLAB 

(total  nodal  points)  *  3  *  NB  (la) 

1 


*  Symbols  used  in  the  WESLAYER  program  are  defined  in  the  input 
guide  (Table  1). 
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or 


(total  nodal  points  in  the  slabs 


(lb) 


where  NB  is  the  half  bandwidth  and  is  equal  to  (NY  +  2)  x  2  .  The 
symbol  NY  represents  the  number  of  nodes  in  the  Y-direction.  When  the 
dimensions  of  C  and  G  are  changed  in  the  main  program,  they  should 
also  be  changed  in  the  corresponding  subroutines. 

16.  The  dimensions  of  many  important  variables  in  the  present 
program  are  limited  and  should  be  increased  if  the  number  of  elements  is 
increased.  The  limited  dimensions  of  variables  are  presented  as  follows: 

a.  The  maximum  nodal  number  (N023)  is  TO  and  the  maximum 
number  of  elements  is  60.  To  increase  element  size,  the 
dimensions  of  the  variables  DF  ,  PPF  ,  PF,  CURL  , 

FO  ,  AB  ,  NCC  ,  REA  ,  and  REACT  should  be  increased 
in  accordance  with  the  total  number  of  nodes  N023  ;  and 
the  dimension  of  the  variable  DEF  is  to  be  increased 
as  the  total  number  of  elements. 

b.  The  maximum  numbers  of  NMCF  and  NPRINT  are  35  and 
30,  respectively.  When  the  value  of  NPRINT  is  changed, 
the  dimension  of  the  variable  NP  should  be  changed 
accordingly. 

c_.  The  dimensions  of  H  (so  is  A),  F  ,  and  CO  are 
determined  by  the  formulas 

N023  x  (N023  +  l)/2;  3  *  N023;  and  3  *  NB  x  NY 

respectively.  Variable  A  is  in  the  subroutine  SMIN  . 

d.  The  maximum  number  of  nodes  that  have  gaps  under  the 

pavement  is  30.  When  this  is  increased,  the  dimensions 
of  the  variables  NG  and  GAP  should  also  be  increased. 

e_.  The  maximum  number  of  loads  is  controlled  by  the  variable 
NL  .  When  the  dimension  of  NL  is  increased,  the  di¬ 
mensions  of  XDA  and  YDA  should  also  be  increased 
accordingly. 

_f.  The  maximum  number  of  layers  is  5*  If  more  layers  are 
desired,  the  dimensions  of  many  variables  in  the  common 
block  RMCOY  should  be  changed. 

Element  and  node  numbering  system 

17.  Beginning  from  the  left  slab  and  ending  at  the  right  slab. 
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the  nodes  and  elements  are  numbered  consecutively  from  bottom  to  top 
and  then  from  left  to  right. 

Symmetries 

18.  The  application  of  the  finite  element  method  for  analyzing 
rigid  pavements  involves  solving  a  large  set  of  simultaneous  equations. 
However,  because  of  symmetry,  the  number  of  simultaneous  equations  could 
be  greatly  reduced  by  considering  only  one-quarter  or  one-half  of  the 
slab.  The  symmetry  is  with  respect  to  the  load,  the  pavement  geometry 
and  property,  the  finite  element  grid  layout,  and  the  load  transfer 
device  along  the  joint.  Users  are  strongly  urged  to  take  advantage  of 
the  symmetry  option  provided  by  the  program  to  arrange  the  loadings  in 
such  a  way  that  the  problem  becomes  symmetrical.  Coded  data  input  for 
symmetrical  example  problems  are  presented  in  Part  V.  It  should  be 
pointed  out  that  symmetry  should  not  be  placed  at  a  joint. 

19.  When  the  effects  due  to  temperature  and  loadings  are  con¬ 
sidered  separately,  the  computed  results  due  to  temperature  alone  are 
expected  to  be  symmetrical  with  respect  to  the  pavement  geometry.  For 
instance,  the  stresses  and  deflections  are  the  same  at  the  four  corner 
nodes  in  a  square  slab  subjected  to  a  temperature  warping.  This  may 
not  be  the  case,  however,  if  the  finite  element  grid  lines  are  not 
divided  symmetrically.  In  practical  cases,  smaller  elements  can  be 
used  around  the  applied  loads,  which  may  result  in  a  nonsymmetrical 
finite  element  grid  pattern.  If  this  is  the  case,  the  computed  results 
due  to  temperature  alone  may  not  be  symmetrical  as  they  should  be  and 
consequently  may  affect  to  a  certain  extent  the  final  results  when  the 
temperature  effect  is  combined  with  that  of  the  load.  The  error  in 
most  cases  is  insignificant  because  the  load  effect  usually  outshadows 
the  temperature  effect.  Nevertheless,  users  should  be  aware  of  this 
possible  discrepancy. 

20.  In  Figure  h  of  Report  2  of  this  series,  the  loads  are  placed 
at  the  pavement's  center  next  to  the  joint.  Smaller  elements  are  used 
around  the  loads  and  larger  elements  are  used  elsewhere.  Although  the 
finite  element  pattern  is  symmetrical  in  the  up-and-down  direction  with 
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respect  to  the  pavement's  center  line  and  symmetrical  in  the  left-and- 
right  direction  of  the  two-slab  pavement  with  respect  to  the  joint,  the 
size  of  each  element  is  not  identical.  Consequently,  if  there  is  no 
moment  transfer  along  the  joint,  the  computed  results  due  to  the  tempera¬ 
ture  effect  at  nodes  1  and  57  are  not  equal,  as  they  theoretically  should 
be.  Consequently,  the  final  computed  results  are  not  strictly  correct. 
However,  the  error  is  believed  to  be  insignificant  when  the  effect  of 
applied  loads  is  combined. 

Half  bandwidth 

21.  The  definition  of  the  half  bandwidth  of  a  matrix  can  be 
found  in  any  structures  book.  The  size  of  the  half  bandwidth  directly 
influences  the  size  of  the  storage  space.  A  proper  nodal  numbering  sys¬ 
tem  may  reduce  the  size  of  the  half  bandwidth.  This  is  illustrated  in 
the  two  different  numbering  systems  shown  in  Figure  6  of  Report  2  of 
this  series.  Both  slabs  in  Figures  6a  and  6b  have  20  nodes  and  12  ele¬ 
ments,  but  the  half  bandwidth  for  the  arrangement  shown  in  Figure  6a  is 
(It  +  2)  «  3  =  18)  and  that  of  Figure  6b  is  (5  +  2)  *  3  =  21.  The  rule 
of  thumb  is  to  arrange  the  finite  element  grid  with  the  side  having 
fewer  nodes  in  the  vertical  direction. 

Weight  of  the  concrete  slab 

22.  In  the  classical  Westergaard  solution,  the  weight  of  the 
slab  is  not  considered  in  the  computation.  The  consideration  of  the 
weight  of  the  slab  is  an  option  in  this  computer  program.  When  tempera¬ 
ture  and  loads  are  not  considered  and  the  subgrade  is  uniform  and  in 
full  contact  with  the  slab,  the  weight  of  the  slab  only  causes  the  slab 
to  settle  uniformly  and  induces  no  bending  in  the  slab.  Consequently, 
stresses  are  not  induced  in  the  slab.  In  some  cases,  the  consideration 
of  the  weight  of  the  slab  is  mandatory,  as  discussed  below. 

23.  The  major  difference  in  procedure  between  full  and  partial 
contact  between  the  slab  and  the  subgrade  is  that  it  is  not  necessary  to 
consider  the  weight  of  the  slab  in  the  case  of  full  contact,  but  the 
weight  of  the  slab  must  be  considered  in  the  case  of  partial  contact. 

2b.  When  problems  involve  only  temperature  warping  (no  exter¬ 
nally  applied  forces),  the  weight  of  the  slab  must  be  considered  to 
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avoid  the  possible  divergence  of  the  solution.  This  is  particularly 
true  when  gaps  exist  under  some  of  the  nodes.  For  the  case  of  partial 
contact,  the  weight  of  the  slab  must  be  considered  even  when  temperature 
is  not  considered. 

Selected  points  of  stress  computations 

25.  While  the  displacements  are  computed  automatically  for  every 
nodal  point,  the  stresses  are  computed  only  on  request.  The  stress 
matrix  is  used  each  time  when  the  stresses  at  a  nodal  point  are  com¬ 
puted.  Some  computer  time  can  be  saved  if  the  stresses  at  only  a  few 
selected  nodes  are  computed. 

Temperature  considerations 

26.  When  temperature  is  considered,  the  dimensions  of  the  two 
slabs  must  be  identical. 

27.  The  computed  initial  curlings  are  independent  of  the  arrange¬ 
ment  of  the  finite  element  grid  pattern.  The  amount  of  initial  curling 
at  each  node  is  computed  by  Equation  10b  in  Report  1  of  this  series. 

The  only  variable  in  Equation  10b  is  the  distance  R  between  the  center 

of  each  slab  to  the  node  where  the  curling  is  computed. 

Correctness  and  divergence 
of  the  obtained  solution 

28.  Users  of  the  computer  program  should  always  be  scrupulous 
with  the  results  computed.  The  stresses  and  deflections  could  be  com¬ 
puted  and  tabulated,  but  the  values  may  not  be  meaningful.  Certain 
features  in  the  program  deserve  special  attention  and  are  explained  in 
the  following  paragraphs. 

29.  Number  of  iterations.  When  the  number  of  iteration  IC 
has  reached  the  maximum  allowable  number  of  iteration  ICL  and  ICLF 
(item  10  of  Table  l),  the  solution  has  not  converged.  The  problem 
should  be  recomputed  with  larger  values  of  ICLF  (and  also  ICL  in 
certain  cases).  However,  it  may  be  wise  at  this  stage  to  see  whether 
the  solution  obtained  is  good  enough  for  engineering  purposes.  In  some 
cases,  a  solution  may  not  be  obtainable  if  the  convergence  criterion  is 
too  strict.  The  same  reasoning  can  be  used  to  check  the  number  of 
iteration  NIC  for  subgrade  contact  against  the  maximum  allowable 
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number  of  iteration  NCYCLE  .  The  value  of  ICL  is  not  as  critical  as 
the  value  of  ICLF  ;  however,  a  large  difference  between  the  actual 
value  of  IC  (printed  in  the  output)  and  the  specified  ICL  is  not 
recommended. 

30.  Reduction  of  relaxation  factor  RFJ  .  If  convergent  re¬ 
sults  cannot  be  obtained,  the  program  reduces  the  factor  automatically. 

A  value  for  the  relaxation  factor  that  is  too  small  results  in  a  shear 
transfer  that  is  too  small  across  the  joint  during  each  iteration;  con¬ 
sequently,  the  computed  results  could  be  erroneous  because  the  con¬ 
vergence  of  the  solution  is  artificially  enforced. 

31.  Symmetries ,  When  symmetry  in  a  given  direct  ion  is  use  i  and 
the  deflections  and  stresses  across  a  certain  Joint  are  supposed  to  be 
equal,  the  efficiency  of  load  transfer  across  the  joint  should  be  input 
only  as  100  percent.  Otherwise,  erroneous  results  will  be  computed. 

32.  Computed  stresses  and  deflections.  The  solutions  obtained 
from  the  finite  element  application  are  by  no  means  completely  correct; 
they  are  merely  close,  acceptable  approximations.  The  correctness  of 
the  computed  larger  values  is  more  important.  The  smaller  values  com¬ 
puted  at  insignificant  locations  of  the  pavement,  such  as  at  placed  far 
away  from  the  load,  are  of  no  significance  in  engineering  problems. 

Input  Guide 

33.  The  input  guide  for  the  program  is  given  in  Table  1,  with 
detailed  explanations  of  each  entry  presented  as  follows: 

a.  Item  1:  Number  of  Runs  Card  (15). 

Notes  Columns  Variables  _ Entry _ 

(l)  1-5  NR UN  Number  of  runs  to  be 

computed 

NOTES: 

(l)  The  number  of  runs  is  first  specified  at  the  onset  of 
computations.  The  nature  of  the  problems  in  each  individual 
run  is  generally  different.  However,  results  of  one  run  can 
be  used  in  the  next  run  immediately  followed  by  the  input 
NREAD  or  NREAD1  .  They  are  explained  in  Item  6. 
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Table  1.  Input  Guide  for  WESLAYER — Slabs  on  Layered  Elastic  Solids 


(Continued)  (Sheet  1  of  5) 


GENERAL  purpose  data  form 


(Continued)  (Sheet  2  of  5) 


GENERAL  PURPOSE  DATA  FORM 


(Continued)  (Sheet  3  of  5) 


Table  1  (Continued) 


(Continued)  (Sheet  It  of  5) 


GENERAL  PURPOSE  DATA  FORM 


Identification  Card  (A8o). 

Notes 

Columns 

Variables 

Entry 

(1) 

1-80 

T1TLE{12a6) 

Enter  the  heading  in¬ 
formation  to  be  printed 
with  the  output 

NOTES: 

(l)  Begin  each  new 

problem  with  a 

new  heading  card. 

Item  3: 

Data  Card  (l6l5). 

Notes 

Columns 

Variables 

Entry 

(1) 

1-5 

NSLAB 

Number  of  slabs  in  the 
model 

6-10 

PR 

Poisson's  ratio  of  the 
concrete,  generally  use 
0.15 

11-15 

T 

Thickness  of  the 
concrete 

16-25 

YM 

Young's  modulus  of  the 
concrete 

(2) 

26-35 

YMS 

Young's  modulus  of  the 
subgrade 

(2) 

36-40 

PRS 

Poisson’s  ratio  of  the 
subgrade 

(3) 

41-45 

NSYM 

Condition  of  symmetry 

(4) 

46-50 

N0TC0N 

Total  number  of  nodes  at 
which  reactive  pressure 
is  initially  set  at  zero 

(5) 

51-55 

NSTORE 

Options  for  thermal 

stress  and  thermal 
deflections : 


EQ.O  The  values  deter¬ 
mined  from  the 
previous  problem 
are  not  used,  or 
it  is  a  new 
problem 

EQ.l  The  values  deter¬ 
mined  from  the 
previous  problem 
are  used 
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Notes 


Columns 


Variables 


56-60  N PUNCH 

(6)  61-65  NB 

66-70  LTR 


_ Entry _ 

Option  for  punching 
values  of  thermal 
stresses  and  deflections 
on  cards: 

EQ.O  No 

EQ.l  Yes 

Half  bandwidth  of  the 
matrix 

Option  for  computing 
load  transfer  along  the 
joint : 

EQ.O  Efficiency  of 

shear  transfer  is 
read  in 

EQ.l  Dowel  information 
is  read  in 


NOTES : 

(1)  The  maximum  number  of  slabs  is  2. 

(2)  When  the  number  of  NLA  is  1,  i.e.,  the  subgrade  is 
modeled  as  a  linear  elastic  solid,  YMS  and  PRS  can  be 
any  number. 


(3)  Assign  1  when  no  symmetry  exists,  2  when  symmetric  with 
respect  to  the  Y-axis,  3  when  symmetric  with  respect  to  the 
X-axis,  4  when  symmetric  with  respect  to  both  the  X-  and 
Y-axes,  and  5  for  four  slabs  symmatrically  loaded.  When 
subgrade  stresses  and  deflections  are  computed,  symmetry 
should  be  used  with  caution.  When  either  NSX  or  NSY 
is  not  zero,  the  total  number  of  nodal  reactive  forces  is 
reduced  to  one-half,  and  when  both  NSX  and  NSY  are  not 
zero,  the  total  number  of  nodal  reactive  forces  is  reduced 
to  one-quarter.  Users  are  urged  to  distinguish  between  a 
four-slab  pavement  symmetrically  loaded  (assign  5  for  NSYM 
for  this  case)  and  a  single  slab  symmetrically  loaded  (as¬ 
sign  4  for  this  case).  In  the  former  case,  the  slabs  are 
divided  by  joints  and  the  shear  transfer  across  the  joints 
can  be  assumed  either  as  100  percent  or  as  by  dowel  bars. 

The  assumption  of  shear  transfer  less  than  100  percent  is 
meaningless  because  of  the  symmetrical  loading.  The  sub¬ 
grade  reactive  force  at  a  node  along  a  neutral  axis  is  only 
one-half  of  the  nodal  force  at  the  node  when  symmetry  is 
not  used.  Symmetry  should  not  be  used  at  nodes  along  a 
joint . 


(4)  If  the  subgrade  soil  at  certain  nodal  points  is  known 
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to  be  not  in  contact  with  the  pavement  due  to  pumping  or 
plastic  deformation,  the  subgrade  reactive  pressure  at  these 
nodes  can  be  initially  set  at  0  to  obtain  speeding  conver¬ 
gence.  If  NCYCLE  =  1  (NCYCLE  is  listed  in  Ttem  It),  these 
nodes  will  never  be  in  contact.  If  NCYCLE  >  1  ,  these 
nodes  may  or  may  not  be  in  contact,  depending  on  calculated 
results . 

(5)  In  pavement  design,  engineers  are  interested  in  stresses 
induced  by  the  applied  load  and  the  temperature  warping.  In 
pavement  research,  however,  engineers  tend  to  measure  only 
stresses  due  to  the  applied  load  because  thermal  stresses 
are  difficult  to  measure.  To  compute  stresses  and  deflec¬ 
tions  by  the  load  alone,  two  separate  but  consecutive  runs 
have  to  be  conducted.  The  first  run  computes  the  thermal 
stresses  alone.  This  is  done  by  setting  NSTORE  =  0  , 

NWT  =  1  ,  NTEMP  =  1  ,  NGAP  >1  (if  this  is  the  case), 
NOTCON  >0  (if  this  is  the  case),  and  NLOAD  =0  in  the 
first  run.  In  the  second  run,  the  stresses  induced  by  the 
applied  load  and  the  temperature  warping  are  computed  by 
setting  NSTORE  =  1  ,  NWT  =  1  ,  NTEMP  =  1  ,  NGAP  >1  (if 
this  is  the  case),  NOTCON  >0  (if  this  is  the  case),  and 
NLOAD  equal  to  the  actual  number  of  loads.  The  differences 
between  those  values  computed  in  the  first  and  second  runs 
are  the  stresses  and  deflections  due  to  the  applied  load 
along.  Note  that  when  temperature  is  considered,  the  slab 
and  the  subgrade  may  be  in  partial  contact,  the  principle 
of  superposition  may  no  longer  be  held  true  (see  paragraph 
47,  Report  1  of  this  series).  It  should  also  be  pointed 
out  that  in  the  case  of  the  first  and  second  runs  discussed 
above,  the  measured  gaps  that  are  input  should  not  include 
the  gaps  due  to  the  temperature  warping  because  they  are 
to  be  computed. 

(6)  The  half  bandwidth  NB  should  be  equal  to  or  greater 
than  (NY  +  2)  *  3  ,  where  NY  is  the  number  of  nodes  in 
the  Y-direction.  It  was  found  that  when  the  number  of  equa¬ 
tions  is  large  while  the  half  bandwidth  is  small,  the  dis¬ 
placements  may  not  converge,  and  a  larger  bandwidth  should 
be  used.  For  NSLAB  =  1  ,  the  solution  can  be  obtained 
without  iteration  if  a  full  bandwidth  NB  is  used. 


Item  4 : 

Data  Card 

(615). 

Notes 

Columns 

Variables 

Entry 

1-5 

NX1 

Number 

X-axis 

of 

in 

nodes  along 
slab  1 

the 

(1) 

6-10 

NX2 

Number 

X-axis 

of 

in 

nodes  along 
slab  2 

the 

11-15 

NY 

Number 

Y-axis 

of 

nodes  along 

the 

Notes 

Columns 

Variables 

Entry 

16-20 

NCYCLE 

Maximum  number  of  cycles 
for  checking  subgrade 
contact,  generally  use 

10  or  more 

(2) 

21-25 

NPRINT 

Number  of  nodes  at  which 
stresses  and  deflections 
are  to  be  printed 

26-30 

NLA 

Number  of  layers  in  the 
subgrade 

31-35 

NCOMP 

OPTIONS  for  computing 

stresses  and  deflections 
in  the  subgrade: 

E^.O  No 

EQ.l  Yes 

NOTES: 

(1)  NX2  =  0  if  NSLAB  =  1  . 

(2)  The  deflections  at  each  node  are  computed  in  the  pro¬ 
gram  but  the  stresses  at  any  node  are  computed  only  on 
request . 

e •  Item  5'-  Stresses  Print  Card  (l6l5).^ 1  ^ 


Notes 

Columns 

Variables 

Entry 

(2) 

1-5 

NP(  I ) 

Nodal  number  whose 
stresses  are  to  be 
printed 

NOTES: 

(l)  Deflections  are  printed  for 

all  nodal  points. 

(2)  Continue  the  input  until  the  number  of  NPRINT 
is  satisfied.  Continue  to  next  data  card  if  NPRINT 
greater  than  l6. 

( Item  ^ ) 
is 

Item  6: 

Sub grade 

Material  Properties  Card  (8F10.2). 

Note: 

Skip  Item  6 

is  NLA  =  1  . 

Notes 

Columns 

Variables 

Entry 

1-10 

E(l) 

Young's  modulus  in 
1 

layer 

11-20 

V(l) 

Poisson's  ratio  in 
1 

1  ayer 

21-30 

E(  2) 

Young's  modulus  in 
layer  2 

25 


Notes 

Columns 

Variables 

Entry 

33  -HO 

V(2) 

Poisson's  ratio  in 

layer  2 

E(NLA) 

V  ( NLA ) 

Item  7 

Subgrade  Layer  Thickness 

Card  (8F10.2). 

Note  : 

Skip  Item  7  if 

NLA  =  1  . 

Notes 

Columns 

Variables 

Entry 

1-10 

HI(  1 ) 

Thickness  of  the  first 
layer 

11-20 

HI  ( 2 ) 

Thickness  of  the  second 
layer 

HI(NN) 

Thickness  of  the  last 
layer,  where  NN  =  NLA 
-  1 

Item  8 

Nodal  Points 

Coordinates 

Cards  (8F10.5). 

Card  1 

:  X-coordinates 

Notes 

Columns 

Variables 

Entry 

(1) 

1-10 

X(l,l) 

X-coordinate  of  the  first 
node  of  slab  1 

11-20 

X(  1 ,2 ) 

X-coordinate  of  the 
second  node  of  slab  1 

X( 1 ,NX1 ) 

X-coordinate  of  the  NX1 
node  of  slab  1 

(2) 

X(2,l) 

X-coordinate  of  the  first 
node  of  slab  2 

X(2,2) 

X-coordinate  of  the 
second  node  of  slab  .? 

X(2,NX2) 

X-coordinate  of  the  NX2 
node  of  slab  2 

Card  2 

:  Y-coordinates 

Notes 

Columns 

Variables 

Entry 

1-10 

Y(l) 

Y-coordinate  of  the  first 

node 
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Notes 


Columns  Variables 

11-20  Y(2) 


_ Entri 

Y-coordinate  of  the 
second  node 


Y(NY) 


Y-coordinate  of  the  NY 
node 


NOTES: 


(1)  The  nodes  at  both  sides  of  the  joint  should  be  identi¬ 
cal;  i.e.,  if  the  slab  length  is  180  in.  and  12  nodes  are 
used  in  the  X-direction,  the  coordinates  of  X(y,12)  and 
X(2,l)  are  both  l80. 

(2)  Skip  input  for  slab  2  if  NSLAB  =  1  . 

When  temperature  is  considered,  the  dimensions  of  two  slabs 
have  to  be  identified.  Otherwise,  erroneous  results  will 
be  computed. 

i.  Item  9:  Subgrade  Contact  Card  (l6l5). 

Note:  Skip  Item  9  if  NOTCON  =  0  ;  i.e.,  the  slabs  are 
initially  in  full  contact  with  the  subgrade. 

Notes  Columns  Variables  _ _ Entry _ 

(l)  NZ(l)  Nodal  number  at  which  re¬ 

active  pressure  is  ini¬ 
tially  assumed  zero 

NOTES : 


(l)  Continue  the  input  until  the  number  of  NOTCON  is 
satisfied.  Continue  to  next  data  card  if  NOTCON  is 
greater  than  l6. 

j_.  Item  10:  Data  Card  (615). 


Notes  Columns 
1-5 


6-10 


11-15 


Variables 

NGAP 


NTEMP 


NLOAD 


_ Entry _ 

Total  number  of  nodes  at 
which  a  gap  exists  be¬ 
tween  slab  and  subgrade ; 
assign  0  if  no  gap  exists 

Condition  of  temperature 
warping : 

EQ.O  Temperature  gra¬ 
dient  is  zero 

EQ.l  Temperature  gra¬ 
dient  is  not  zero 

Number  of  elements  on 
which  load  is  applied; 
use  0  if  there  is  no  load 
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Columns 

16-20 


Notes 

(1) 


(2) 


(2) 


(3) 


(10 


21-25 


26-30 


31-35 


36-U0 


U1-U5 


Variables 

NMCF 


ICL 


TCLF 


NCR 


NWT 


IGNOR 


_ Entry _ 

Number  of  concentrated 
nodal  forces  and  moments 
that  are  to  be  read  in; 
assign  0  if  no  moments  or 
forces  are  applied 

Maximum  number  of  itera¬ 
tions  allowed  for  coarse 
control;  generally  use 
U9  or  more 

Maximum  number  of  itera¬ 
tions  allowed  for  fine 
control;  generally  use 
199  or  more 

Total  number  of  nodal 
points  for  checking 
convergence 

Weight  of  slab  considera¬ 
tion  : 

EQ.O  Weight  is  not 
considered 

EQ.l  Weight  is  con¬ 
sidered 

A  parameter  indicating 
whether  the  reduction  of 
relaxation  factor  RFJ 
should  be  ignored: 


EQ. 

,0 

If 

RFJ  is  reduced 

EQ. 

,1 

If 

RFJ  is  not  re- 

duced  whenever  the 
results  diverge 


NOTES : 

(l)  The  concentrated  force  is  considered  to  be  positive  if 
it  is  acting  downward  and  to  be  negative  if  it  is  acting 
upward.  Positive  moment  follows  the  right-hand  screw  sys¬ 
tem,  as  shown  in  Figure  1  of  Report  1  of  this  series.  The 
program  is  dimensioned  for  50  concentrated  forces  and  mo¬ 
ments.  If  NMCF  is  greater  than  50,  dimensions  of  NFF  , 
NFI  ,  and  NF  must  be  increased. 


(2)  See  note  1  of  Item  11.  Coarse  and  fine  controls  are 
used  before  and  after  the  subgrade  contact  condition  is 
determined.  For  a  given  contact  condition,  coarse  control 
is  used  to  ch.ck  the  deflection  criterion.  Once  the  sub¬ 
grade  contact  condition  is  finally  determined,  fine  control 
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is  used  to  obtain  accurate  solutions.  In  the  case  where 
NCYCLE  =  1  ,  coarse  control  is  still  used  prior  to  the  use 
of  fine  control.  Note  ICLF  should  always  be  greater  than 
ICL  . 

In  some  problems,  ICLF  may  be  exhausted  before  the  cri¬ 
terion  DELF  is  satisfied.  Before  rejecting  the  solution, 
it  may  be  wise  to  check  to  see  how  far  the  solution  is  from 
satisfying  the  criterion.  For  instance,  if  DELF  =  0.001 
and  computed  convergence  is  0.002  or  0.0025  and  the  computed 
results  seem  to  be  reasonable,  the  solution  may  be  con¬ 
sidered  acceptable.  In  some  problems,  it  may  be  very  hard 
to  satisfy  the  specified  convergence  criterion. 

(3)  Two  or  three  points  should  suffice.  They  should  be 
selected  under  or  near  the  load  where  stresses  and  deflec¬ 
tions  are  the  greatest.  When  NSLAB  =  2  ,  the  convergence 
points  should  be  selected  along  the  joint. 

(L)  The  use  of  ICNOR  is  to  increase  the  flexibility  of 
the  program.  In  some  cases,  it  may  be  desirable  to  check 
the  convergence  condition  when  the  relaxation  factor  is 
fixed  at  a  certain  value.  It  was  found  that  in  some  cases 
the  computed  results  are  better  when  the  problem  is  com¬ 
puted  with  a  constant  RFJ  of  0.5  than  when  computed  with 
either  a  smaller  constant  RFJ  or  a  varying  RFJ  (reducing 
automatically).  For  instance,  if  the  loads  are  symmetri¬ 
cally  placed  on  the  slab,  the  stresses  should  be  symmetrical 
at  corresponding  points.  It  was  found  that  the  magnitudes 
of  the  stresses  computed  at  the  corresponding  points  were 
much  closer  when  computed  with  a  constant  RFJ  of  0.5  than 
when  computed  otherwise.  It  is  thus  recommended  to  use 
RFJ  =  0.5  or  0.1  and  IGNOR  =1  in  the  computation  if  the 
solution  can  be  obtained.  The  use  of  an  RFJ  that  is  too 
small  can  artificially  accelerate  the  convergence,  but  the 
results  obtained  may  be  erroneous. 

k.  Item  11:  Data  Card  (8F10.R). 

Notes  Columns  Variables  _ Entry _ 

1-10  TEMP  Difference  in  tempera¬ 

ture,  in  degrees  Fahren¬ 
heit,*  between  top  and 
bottom  of  slab: 

EQ. positive  slab  curled 
upward 

EQ. negative  slab  curled 
downward 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (ST)  units  is  presented  on  page  3. 
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Notes 

Columns 

Variables 

Entry 

11-20 

Q 

Uniformly  applied  pres¬ 
sure  in  psi 

(1) 

21-30 

DEL 

Tolerance  of  convergence 
for  coarse  control;  usu¬ 
ally  use  0.01 

(1) 

31-^40 

DELF 

Tolerance  of  convergence 
for  fine  control ;  usu¬ 
ally  use  0.001 

(2) 

141-50 

RFJ 

Initial  relaxation  factor 
at  the  joint ;  generally 
use  0.5  or  0.1 

(3) 

51-60 

EFF 

Efficiency  of  shear 
transfer  across  the 
joint.  It  ranges  from 

1  to  0 

NOTES : 

(1)  DEL 

and  DELF 

correspond  to 

ICL  and  ICLF  in 

Item  10,  respectively. 

(2)  The  use  of  RFJ  is  to  reduce  the  deflections  trans¬ 
ferred  from  one  slab  to  the  other.  The  program  reduces  the 
value  of  RFJ  automatically  when  the  solution  diverges. 

When  dowel  bars  are  used,  it  was  found  that  the  solution 
can  be  obtained  only  when  RFJ  is  reduced  to  0.01.  When 
dowel  bars  are  not  used,  however,  an  erroneous  solution  was 
obtained  when  RFJ  was  set  to  0.01. 

(3)  The  value  of  efficiency  across  the  joint  varies  from  0 
to  1.  If  LTR  (input  in  Item  3)  is  equal  to  1,  EFF  must 
be  input  as  1.  However,  it  does  not  mean  that  100  percent 
shear  transfer  is  used  in  the  program. 

No  minimum  value  of  RFJ  is  used  in  the  program.  However, 
a  value  too  small  would  result  in  a  shear  transfer  across 
the  joint  that  is  too  small  during  each  iteration,  and  con¬ 
sequently  the  computed  results  could  be  erroneous  because 
the  convergence  of  the  solution  is  "artifically"  enforced. 

When  the  pavement  geometry  and  loading  conditions  in  the 
pavement  are  such  that  the  computed  results  at  corresponding 
locations  in  two  slabs  are  supposed  to  be  symmetrical.  Joint 
efficiency  other  than  100  percent  should  not  be  used.  Other¬ 
wise,  erroneous  results  will  be  computed. 

Item  12:  Concentrated  Forces  and  Moments  Card  Qi(I5, 

13,  F10 . 2 )3~ 

Note:  Skip  Item  22  if  NMCF  =  0  . 
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Notes 

Columns 

Variables 

Entry 

1-5 

NF(I) 

Nodal  number  at  which 
concentrated  forces  or 
moments  are  specified 

(1) 

6-10 

NFF(I) 

Nature  of  specified  force 
at  node  I 

(2) 

11-20 

fo[nf(i)-i] 

x  3  +  NFF(I) 

Concentrated  force  or  mo¬ 
ment  at  node  I 

NOTES : 

(l)  NFF(l)  =  1  for  vertical  force,  2  for  moment  about 
X-axis,  and  3  for  moment  about  Y-axis. 

(2)  The  magnitude  of  concentrated  force  or  moment  at  each 
node  input  in  the  equation  number  is  related  to  nodal  number 
I  by  NF(l)  -  1  x  3  +  NFF(l)  .  For  instance,  if  a  moment 
about  Y-axis  is  applied  at  node  13,  the  equation  number  will 
be  ( 13  -  1 )  *  3  +  3  =  39  •  Note  that  the  nodes  are  num¬ 
bered  consecutively  from  bottom  to  top  and  then  from  left 
to  right  beginning  from  the  first  slab  and  ending  at  the 
last  slab. 

Item  13: 

Check  Convergence  Card 

( l6l5 ) . 

Notes 

Columns 

Variables 

Entry 

1-5 

NODCK(l) 

Nodal  number  of  the  first 
node  at  which  the  conver¬ 
gence  is  checked 

6-10 

N0DCK(2) 

Nodal  number  of  the  sec¬ 
ond  node  at  which  the 
convergence  is  checked 

NODC(NCK) 

Nodal  number  of  the  NCK 
node  at  which  the  con¬ 
vergence  is  checked 

Item  lH 

:  Dowel  Bar 

Information 

Card  (5F10.5,  E10.3). 

Note: 

Skip  Item  lb 

if  LTR  =  0 

. 

Notes 

Columns 

Variables 

Entry 

1-10 

BD 

Bar  diameter 

11-20 

BS 

Bar  spacing 

21-30 

WJ 

Joint  opening 

31-U0 

PRSB 

Poisson's  ratio  of  the 
bar 

31 
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Notes 

Columns 

Variables 

Entry 

Ul-50 

YMSB 

Young' s 
bar 

modulus  of  the 

51-60 

DSM 

Modulus 

of  dowel  support 

o_.  Item  15:  Gaps  Read  In  Cards. 

Note:  Skip  Item  15  if  NGAP  =  0  . 
Card  1:  Nodal  Points  (1615) 

Notes  Columns  Variables 
1-5  NG(1) 

6-10  NG(2) 


NG(NGAP) 

Card  2:  Gaps  (8FIO.5) 

Notes  Columns  Variables 
1-10  CURL(NG(l ) ) 

11-20  CURL(NG(2)) 

CURL(NG(NGAP) )  Amount  of  gap  at  node 
NG(NGAP) 

£.  Item  16:  Total  Uniformly  Applied  Load  Card  (F12.2). 

Note:  Skip  Item  16  if  NLOAD  =  0  . 


Notes 

Columns 

Variables 

Entry 

(1) 

1-12 

RLOAD 

Total  uniformly  applied 

load  on  the  slab 

NOTES: 

(l)  The  total  load  refers  to  the  uniformly  applied  load 
only.  The  total  load  should  be  divided  by  2  or  U  if  it  is 
symmetric  with  respect  to  one  axis  (X-  or  Y-axis)  or  both 
the  X-  and  the  Y-axis,  respectively.  Additional  point  loads 
applied  at  nodal  points  are  excluded. 


_ Entry _ 

Nodal  number  of  the  first 
node  at  which  a  gap 
exists 

Nodal  number  of  the  sec¬ 
ond  node  at  which  a  gap 
exists 


Nodal  number  of  the  NGAP 
node  at  which  a  gap 
exists 


_ _ Entry _ 

Amount  of  gap  at  node 
NG(  1 ) 

Amount  of  gap  at  node 
NG(  2 ) 
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c^.  Item  17:  Loading  Cards. 


Note:  Skip  Item  IT  if  NLOAD  =  0  . 
Card  1:  Loading  Elements  (1615; 


Notes 

Columns 

Variables 

Entry 

(1) 

1-10 

NL(1) 

Element  number  of  the 
first  element  over  which 
load  is  applied 

(1) 

11-20 

NL(2) 

Element  number  of  the 
second  element  over  which 
load  is  applied 

NL( NLOAD) 

Element  number  of  the 
last  element  over  which 
load  is  applied 

NOTES: 

(l)  Beginning  from  the  first  slab  and  ending  at  the  last 
slab,  the  nodes  and  elements  are  numbered  consecutively 
from  bottom  to  top  and  then  from  left  to  right . 

Card  2: 

Load  Magnitude  (1+F10.5) 

Notes 

Columns 

Variables 

Entry 

(1) 

1-10 

XDA( I ,1 ) 

Lower  limit  of  loaded 
area  in  element  I  in  X- 
direction 

(1) 

11-20 

XDA ( I , 2 ) 

Upper  limit  of  loaded 
area  in  element  I  in  X- 
direction 

(2) 

21-30 

YDA( I ,1 ) 

Lower  limit  of  loaded 
area  in  element  I  in  Y- 
direction 

(2) 

31-1*0 

YDA( I ,2 ) 

Upper  limit  of  loaded 
area  in  element  I  in  Y- 
direction 

NOTES : 

( 1 )  Use 
element. 

-1  to  +1 

if  the  load  covers  the  whole  length  of  the 

(2)  Use 
element . 

-1  to  +1 

if  the  load  covers  the  whole  width  of  the 

Item  18: 

Subgrade 

Stresses  Card. 

Note:  Skip  Item  18  if  NCONF  =  0  . 
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Card  1: 

Number  of  Computations 

(2110) 

Notes 

Columns 

Variables 

Entry 

1-10 

NZZ 

Number  of  depths  to  be 
computed 

11-20 

NR 

Number  of  offsets  at  each 
depth  to  be  computed 

Card  2: 

Depth  Card 

(8F10.5) 

Notes 

Columns 

Variables 

Entry 

1-10 

ZZ(1) 

Depth  of  first  computa¬ 
tion 

11-20 

ZZ(2) 

Depth  of  second  computa¬ 
tion 

ZZ(NZ) 

Depth  of  the  last  com¬ 
putation 

Card  3: 

Offset  Card 

(8F10.5) 

Notes 

Columns 

Variables 

Entry 

(1) 

•1-10 

XR(1) 

X-coordinate  of  first 
computation 

11-20 

YR(1) 

Y-coordinate  of  first 
computation 

21-30 

XR(2) 

X-coordinate  of  second 
computation 

31-1*0 

YR(2) 

Y-coordinate  of  second 
computation 

XR(NR) 

X-coordinate  of  last 
computation 

YR(NR) 

Y-coordinate  of  last 
computation 

NOTES: 

(l)  Computations  at  each  offset  point  are  made  at  all  the 
high  depths.  The  origin  of  the  coordinates  is  at  nodal 
point  1;  i.e.,  model  of  slab  1.  Refer  to  the  nodal  numbers 
shown  in  Figure  2  of  Report  2  if  stresses  and  deflections 
are  to  be  computed.  The  first  location  is  directly  under 
node  1,  the  second  location  is  midpoint  between  nodes  5 


and  8,  and  the  third  location  is  at  the  center  of  element 
13.  The  input  values  of  XR(l)  ,  YR(l)  ,  XR(2)  ,  YR(2)  , 
XR(3)  ,  and  YR(3)  should  then  be  0,  0,  135,  90,  -135,  and 
-135-  Note  that  nodes  1,  l6,  21,  and  3 6  share  the  same 
location. 


PART  V:  EXAMPLE  PROBLEMS 


3*+.  In  this  Part,  the  input  data  of  four  example  problems  are 
presented.  Printouts  of  the  computer  outputs  for  the  example  problems 
are  also  presented  and  explained. 


Example  Problem  1:  A  Two-Slab  Pavement 
System  on  Elastic  Solids 


35.  Figure  2  shows  the  finite  element  grid  of  a  two-slab 
pavement.  The  nodes  and  elements  are  numbered  consecutively  from  bottom 
to  top  and  then  from  left  to  right.  The  input  data  consist  of  th^  fol¬ 
lowing  information: 

a.  The  concrete  slab  is  8  in.  thick  with  a  Young's  modulus 
of  l*, 000, 000  psi  and  a  Poisson's  ratio  of  0.15.  The 
elastic  subgrade  soil  has  a  Young.' s  modulus  of  10,000 
psi  and  a  Poisson's  ratio  of  0.1+ . 

b.  The  pavement  is  subjected  to  a  uniformly  distributed 
square  load  applied  at  the  center  of  the  slabs.  The 
option  of  symmetry  with  respect  to  the  X-axis  is  used. 

c_.  Dowel  bars  are  not  used  in  the  transverse  joint  con¬ 
necting  the  two  slabs.  The  joint  is  assumed  to  have  a 
certain  percentage  of  efficiency  in  shear  transfer 
across  the  joint. 

d.  The  input  data  for  Example  problem  1  are  given  in 

Table  2.  The  joint  efficiency  in  this  particular  case 
is  100  percent;  i.e.,  the  deflections  at  both  sides  of 
the  joint  are  equal. 

36.  Table  3  shows  the  printout  of  the  computer  output  for  Example 
problem  1.  The  printout  is  in  many  places  self-explanatory.  For  con¬ 
venience  of  explanation,  entry  numbers  are  used  where  explanations  are 
needed. 

Entry  1 

37.  Because  the  option  of  symmetry  with  respect  to  the  X-axis 
was  used,  the  parameter  NSYM  was  input  as  3.  Because  the  option  of 
shear  transfer  across  the  joint  is  used,  the  parameter  LTR  was  input 
as  0.  The  half  bandwidth  NB  =  21  was  the  minimum  required.  In  this 
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Table  2.  Example  Problem  1,  Input  Data  for  a  Two-Slab  Pavement  on  an  Elastic 
Solid,  100  Percent  Shear  Transfer  Across  the  Joint 
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Table  3  (Concluded) 


RFJ  =0.1  is  fixed  as  a  constant 


special  case,  the  relaxation  factor 
by  setting  the  parameter  IGNOR  =  1  . 

Entry  2 

38.  The  initial  curling  and  gap  at  all  nodes  are  0  because  the 
temperature  and  the  gap  are  not  considered  in  the  computation. 

Entry  3 

39-  Since  NCYCLE  =1  ,  i.e.,  the  slab  is  assumed  to  be  in  con¬ 
tact  with  the  subgrade  all  the  time,  NIC  will  not  be  iterated  greater 
than  1. 

Entry  1* 

1*0.  The  deflections  at  each  node  are  tabulated  whether  con¬ 
vergence  requirements  are  met.  Positive  deflection  indicates  downward 
movement . 

Entry  5 

1*1.  The  stress  components  at  preselected  points  are  printed. 
Positive  stress  indicates  that  the  slab  has  compression  at  the  top  and 
tension  at  the  bottom  and  negative  stress  indicates  compression  at  the 
bottom  and  tension  at  the  top.  The  symbols  of  STRESS  XY,  MAJOR, 

MINRO,  and  SHEAR  stand  for  shear  stress,  major  principal  stress,  minor 
principal  stress,  and  maximum  shear  stress,  respectively. 

Entry  6 

b 2.  The  stresses  and  displacements  are  computed  for  one  more 
iteration  for  inspection  of  convergence  by  the  user.  When  the  solution 
correctly  converges,  the  differences  in  the  computed  results  between 
two  iterations  should  be  very  insignificant.  Otherwise,  the  solution 
is  not  convergent. 

1*3.  Because  of  the  existence  of  symmetry  with  respect  to  the 
Y-axis  along  the  transverse  joint,  the  stresses  and  deflections  at  cor¬ 
responding  points  in  the  left  and  right  sides  should  be  equal.  Printed 
results  in  Entry  6  show  that  the  stresses  and  deflections  at  nodal 
points  11,  16,  26,  31,  32,  33,  3l*,  and  35  are  very  close  to  those  com¬ 
puted  in  points  56,  51,  1*1,  36,  37,  38,  39,  and  1*0,  respectively. 

1*1*.  The  following  two  special  features  in  Example  problem  1 
should  be  pointed  out: 


1*2 


su  Because  the  slabs  are  symmetrically  loaded,  the  use  of 
a  shear  transfer  across  the  joint  less  than  100  percent 
will  produce  meaningless  results. 

b.  While  the  system  is  symmetrically  loaded,  the  use  of  the 
symmetry  option  with  respect  to  both  the  X-  and  Y-axis 
is  wrong  because  the  axis  o'*  symmetry  in  the  Y-direction 
lies  on  the  joint  which  is  formidable. 

15.  The  results  shown  in  Table  3  were  computed  with  the  relaxa¬ 
tion  factor  RFJ  =0.1  . 

U 6 .  Table  1  is  another  printout  of  the  computer  output  with  the 
same  conditions  as  those  shown  in  Example  problem  1,  except  that  the 
efficiency  of  shear  transfer  across  the  joint  is  assumed  to  be  80  per¬ 
cent  and  the  loads  applied  at  the  right  slab  are  removed.  Eighty  per¬ 
cent  shear  transfer  denotes  that  the  deflections  at  the  unloaded  (or 
less  heavily  loaded)  slab  are  80  percent  of  the  loaded  (or  more  heavily 
loaded)  slabs.  The  computed  deflections  at  nodes  31  and  36  indicate 
that  the  computed  results  are  correct.  It  should  be  noted  that  since 
80  percent  shear  transfer  is  assumed,  meaningless  results  will  be  com¬ 
puted  if  the  slabs  are  symmetrically  loaded  with  respect  to  the  trans¬ 
verse  joint. 


Example  Problem  2 :  A  Two-Slab  Pavement  System 
on  Elastic  Solids,  Dowel  Bars  Across  the  Joint 

17-  The  same  configuration  of  pavements  used  in  Example  prob- 
men  1  (Figure  2)  is  used  in  this  one.  Dowel  bars  are  used  across  the 
transverse  joint;  the  bars  have  a  diameter  of  1  in.  and  are  placed 
12  in.  apart.  The  input  data  and  the  printout  of  the  computer  output 
are  given  in  Tables  5  and  6,  respectively.  When  dowel  bars  are  used, 
the  options  LTR  and  EFF  are  both  input  as  1.  The  computed  results 
show  that  the  stress  a  and  deflections  in  the  loaded  slab  are 

y 

greater  than  those  in  the  unloaded  slab. 

Example  Problem  3:  A  Two-Slab  Pavement  System 
on  Layered  Elastic  Solids 

18.  Table  7  shows  the  input  data  for  a  two-slab  pavement  system 
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Table  U .  Computer  Output  for  Example  Problem  1,  Elastic  Solid  and 
80  Percent  Shear  Transfer  Across  the  Joint,  One  Loaded  Slab 
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Table  5 •  Example  Problem  2,  Input  Data  for  a  Two-Slab  Pavement  on  an  Elastic 
Solid,  1-in.  Dowel  Bars  Across  the  Joint 
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Table  8 


Printout  of  the  Computer  Output  for  Example  Problem  3, 
Three-Layer  Elastic  Subgrade,  100  Percent  Shear 
Transfer  Across  the  Joint 
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resting  on  a  three-layer  linear  elastic  subgrade.  One  hundred  percent 
joint  efficiency  is  assumed.  The  finite  element  grid  pattern  is  the 
same  as  shown  in  Figure  2. 

49.  Table  8  shows  the  printout  of  the  computer  output  for 
Example  problem  3.  The  layered  subgrade  information  indicates  that  the 
elastic  moduli  of  the  three-layer  subgrade  are  60,000,  30,000,  and  10,000 
psi  and  that  the  Poisson's  ratio  in  each  layer  is  0.4.  Because  of  the 
assumption  of  100  percent  shear  transfer  across  the  joint  and  the  sym¬ 
metrical  loading,  the  stresses  and  deflections  at  corresponding  nodes 
across  the  joint  are  nearly  identical. 

Example  Problem  4 :  A  Four-Slab  Pavement  System  on 
Elastic  Solids,  Symmetrically  Loaded,  with 
100  Percent  Shear  Transfer  Across  the  Joint 


50.  Figure  3  shows  the  finite  element  layout  for  a  four-slab 
pavement  system  with  the  load  applied  symmetrically  at  the  pavement's 
center.  This  is  the  only  option  for  the  analysis  of  a  four-slab 
pavement  system  by  use  of  the  WESLAYER  program.  Although  it  is  a  four- 
slab  system,  the  program  analyzes  only  one-quarter  of  the  system,  i.e., 
one  slab.  Consequently,  100  percent  joint  efficiency  should  be  assumed 


_ I _ 

Figure  3-  Finite  element  layout  for  Example  problem  4, 
symmetrically  loaded  four-slab  system 
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for  the  problem.  The  coordinates  and  loading  condition  input  into  the 
program  are  shown  in  Figure  3. 

51.  The  input  data  for  this  problem  are  given  in  Table  9-  It 
should  be  pointed  out  that  the  parameter  NSYM  is  input  as  5  for  a 
four-slab  system.  The  number  of  slab  NSLAB  is  input  as  1,  instead  of 
U,  because  there  is  only  one  slab  involved  in  the  computation,  (in  the 
case  for  a  single  slab  symmetrically  loaded,  NSYM  should  be  input 

as  4 . ) 

52.  The  printout  of  the  computer  output  for  Example  problem  U  is 
presented  in  Table  10.  The  printout  is  self-explanatory.  Because  of 
symmetry,  the  computed  stresses  and  deflections  at  nodal  points  2,  3,  t, 
and  5  should  be  very  close  to  those  computed  at  nodes  7,  13,  19,  and  25, 
respectively.  The  computed  results  show  that  it  is  true  only  at  nodal 
points  2  and  7-  At  nodes  closer  to  the  free  edge  of  the  pavement,  the 
discrepancy  becomes  greater.  This  is  believed  to  be  caused  by  the  fact 
that  the  slabs  have  different  dimensions  in  the  X-  and  Y -directions .  It 
is  also  to  be  noted  that  at  nodes  1  and  6,  i.e.,  nodes  at  the  corners  of 
the  slabs,  the  stresses  vanish. 


58 


,ple  Problem  U,  Input  Data  for  a  Four-Slab  Pavement  on  an  Elastic 
Symmetrically  Loaded  on  Both  X-  and  Y-Axes,  100  Percent 
Shear  Transfer  Across  the  Joints 


Table  10.  Printout  of  Computer  Output  for  Example  Problem  a  Four-Slab  Pavement 
System  on  Elastic  Solids,  100  Percent  Shear  Transfer  Across  the  Joint 
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Table  10  (Concluded) 


PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


53.  The  computer  program  WESLAYER  has  the  capability  of  ob¬ 
taining  solutions  for  rigid  multicomponent  pavements  with  discontinui¬ 
ties.  The  foundation  soil  can  be  either  elastic  solid  or  layered 
elastic  solid.  Because  of  the  length  of  computer  time  required  and  the 
limitation  of  two  slabs,  it  is  recommended  that  the  program  be  used  for 
research  and  analysis  purposes. 
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